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Abstract 


The  temperature  increase  in  tire  main  body  portion  of  a  kinetic  energy 
penetrator  flying  at  sea  level  with  a  speed  of  1,500  m/s  was  calculated  using  an 
analytical  conduction  analysis  with  a  convective  boundary  condition.  The 
penetrator  was  modeled  as  a  smooth  cylinder  of  uniform  diameter  and  material 
properties.  The  results  from  a  previously  published  computational  fluid 
dynamics  simulation  were  used  to  provide  the  convective  heat  transfer 
coefficient  and  the  temperature  of  the  gas  flowing  over  the  surface  of  the 
penetrator.  A  separation  of  variables  solution  was  used  for  the  parabolic  time- 
dependent  conduction  equation.  The  derived  solution  along  with  material 
properties  for  tungsten  and  depleted  uranium  were  used  to  obtain  the 
temperature  increase  profile  through  the  radius  of  the  penetrator  as  a  function  of 
flight  time.  Both  materials  showed  significant  heating  within  3  s  of  flight.  The 
tungsten  penetrator  exhibited  nearly  uniform  heating  across  the  radius,  whereas 
heating  of  the  depleted  uranium  penetrator  was  confined  to  the  outermost 
region  of  the  cylinder  due  to  its  low  coefficient  of  thermal  conductivity. 
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1.  Introduction 


A  kinetic  energy  (KE)  penetrator  launched  from  a  cannon  travels  to  its  intended 
target  typically  at  speeds  in  excess  of  1,500  m/s.  At  these  velocities,  the  adiabatic 
wall  temperature  of  the  air  reaches  approximately  1300  K  (Guidos  and  Weinacht 
1993),  establishing  the  potential  for  increasing  the  temperature  of  the  penetrator.  If 
the  temperature  of  the  penetrator  is  significantly  increased,  the  mechanical  behavior 
of  the  penetrator  material  will  change,  thereby  affecting  terminal  ballistic 
performance  (Magness  and  Lopatin  1993).  Two  metallic  alloys  commonly  used  for 
KE  penetrators  are  tungsten  heavy  alloys  (WHA)  that  generally  contain  greater  than 
90%  (by  weight)  tungsten  and  depleted  uranium  (DU)  alloyed  with  less  than  1% 
titanium.  The  yield  stress  of  these  materials  has  been  observed  to  decrease 
exponentially  with  temperature  (Zurek  and  Follansbee  1995).  Fracture  properties 
and  fracture  mode  are  temperature  sensitive  as  well.  In  particular,  the  ductile-brittle 
transition  temperature  for  WHA  is  at  room  temperature  or  slightly  above, 
depending  upon  composition,  stress  state,  and  strain  rate  (Larson  and  DeLai  1980; 
Gurwell  1986).  The  ductile-brittle  transition  temperature  of  unalloyed  uranium  is 
0  °C  (Eckelmeyer  1991),  with  the  transition  temperature  of  DU  expected  to  be 
slightly  lower.  An  important  consequence  of  this  transition  temperature  is  that  the 
mode  of  fracture  during  target  interaction  will  change  depending  on  the  penetrator 
flight  distance.  More  significantly,  the  terminal  ballistic  efficacy  of  a  KE  penetrator 
may  depend  on  the  penetrator-target  engagement  distance.  For  these  reasons,  it  is 
necessary  to  know  the  penetrator  temperature  history  during  flight. 

The  flow  field  over  a  KE  penetrator  moving  at  supersonic  and  hypersonic  speeds 
was  determined  by  Guidos  and  Weinacht  (1993)  using  a  three-dimensional, 
parbolized  Navier-Stokes  computational  technique.  Their  numerical  simulations 
provided  perfect  gas  heat  transfer  parameters  for  a  penetrator  of  geometry  similar  to 
the  M829.  Penetrator  speed  was  1,500  m/s.  The  adiabatic  wall  temperatures  and 
the  convective  heat  transfer  coefficients  reported  by  Guidos  and  Weinacht  provide 
the  boundary  conditions  needed  for  transient  conduction  analyses  to  quantify  the 
penetrator  temperature  history.  The  results  contained  within  this  report  were 
obtained  by  solving  a  one-dimensional  transient  heat  conduction  equation  that  uses 
the  Guidos  and  Weinacht  boundary  conditions.  Only  the  cylindrical  main  body 
portion  of  the  penetrator  was  considered  for  the  conduction  solution. 
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2.  Conduction  Model  and  Solution  Technique 


2.1  One-Dimensional  Conduction  Model 

The  heat  transfer  model  used  to  describe  the  temperature  increase  of  a  KE 
penetrator  traveling  through  shock-heated  air  consists  of  a  governing  equation 
derived  from  the  first  law  of  thermodynamics  and  Fourier's  law,  and  a  set  of  initial 
and  boundary  conditions.  The  general  conduction  model  and  its  solution  can  be 
found  in  advanced  heat  conduction  texts  (see  Poulikakos  [1994]  for  example).  It  is 
derived  here  for  completeness  and  to  fully  document  the  assumptions  and  modeling 
simplifications  that  were  made  and  the  solution  technique  used.  The  heat 
conduction  equation  obtained  from  the  rate  form  of  energy  conservation  is  given  as 

^-av’r+e,.,  (1) 


where  T  is  temperature,  t  is  time,  Qgen  is  the  rate  of  internal  heat  generation  per  unit 
volume,  and  a  is  the  thermal  diffusivity  defined  as  a  =  k/pc(k  is  the  thermal 
conductivity,  p  is  the  bulk  density,  and  c  is  the  specific  heat).  Note  that 
compressibility  effects  are  neglected  in  equation  1,  implying  p  =  constant  and 
cp  =  cv  =  c  (where  cp  and  cv  are  the  specific  heat  at  constant  pressure  and  volume, 

respectively),  and  that  the  thermal  conductivity  of  the  material  is  assumed  constant. 
As  will  be  seen,  it  shall  be  convenient  to  work  in  a  cylindrical  coordinate  system 
(r,  6,  z).  In  cylindrical  coordinates,  equation  (2)  has  the  following  form: 


1 8T  a2r  i  d2r  a2rl  • 

r  dr  dr 2  r2  dd2  dz 2  +  sen 


(2) 


Equation  2  can  be  simplified  by  considering  the  geometry  of  a  typical  KE  penetrator, 
shown  in  Figure  1.  The  penetrator  main  body  was  simplified  to  a  right-circular 
cylinder  of  length  L  and  radius  R,  with  a  smooth  surface  (as  was  also  done  in  the 
computational  model  of  Guidos  and  Weinacht  [1993]).  Since  the  focus  of  this  effort 
is  to  identify  the  temperature  history  of  the  main  body,  effects  from  the  nose  cap  and 
fin  were  neglected.  This  leads  to  the  penetrator  configuration  and  heating  process 
depicted  in  Figure  2,  which  is  amenable  to  description  in  a  cylindrical  coordinate 
system.  Noting  symmetry  in  the  0  plane,  the  transient  heat  conduction  is  at  most 
two-dimensional,  occurring  in  the  r-z  plane.  Scaling  arguments  show  that 
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Figure  1.  Geometry  of  a  typical  KE  penetrator. 


Figure  2.  Simplified  penetrator  for  heat  transfer  model. 
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Since  L»R,  equations  3  and  4  imply  that 
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(5) 


Therefore,  the  heat  conduction  equation  expressed  in  cylindrical  coordinates  (2) 
reduces  to  the  following  one-dimensional  governing  equation  (neglecting  internal 
heat  generation): 

1  dT  1  8T  d2T 

- = - +  — r-.  (6) 

a  dt  r  dr  dr 


One  initial  condition  and  two  boundary  conditions  are  required  for  solving 
equation  6.  The  initial  condition  is  a  statement  that  the  penetrator  is  at  ambient 
temperature  when  f  =  0.  The  first  boundary  condition  is  recognizing  that  there  is  no 


heat  flux  at  the  centerline  of  the  penetrator,  i.e.,  a  symmetry  condition.  The  final 
boundary  condition  is  obtained  by  noting  that  if  a  body  is  in  contact  with  a  moving 
convective  fluid,  the  heat  conduction  in  the  body  equals  the  heat  convected  by  the 
fluid.  Mathematically,  the  initial  and  boundary  statements  are  the  following: 

f  =  0 :  T  =  To,  (7) 


r=0:^  =  0, 

dr 


and 


r  =  R:~k—  =  h(r-TJ,  (9) 

where  To  is  the  initial  temperature  of  the  penetrator,  Tw  is  the  temperature  of  the 
shock  heated  air,  and  h  is  the  convective  heat  transfer  coefficient.  The  method  of 
separation  of  variables  will  be  used  to  solve  the  one  dimensional  conduction 
equation  6.  For  this  method  to  apply,  the  governing  equation  as  well  as  both 
boundary  conditions  must  be  homogeneous.  While  equations  6  and  8  are 
homogeneous,  equation  9  is  not.  To  remedy  this  situation,  the  following  variable 
substitution  is  used: 


5© 

r  =  R:-k - =  h&.  (14) 

dr  v  ’ 

Proceeding  with  the  separation  of  variables  solution,  let  ©  be  of  the  product  form 

e(r,t)  =  F(r)G(t),  (15) 


where  the  functions  F(r )  and  G(t)  are  to  be  determined.  Putting  the  assumed  form  of 
©  from  equation  15  into  the  governing  equation  11  and  performing  the  standard 
separation  of  variables  analysis  (see  Greenberg  [1988]  for  additional  details)  yields 
the  following  two  ordinary  differential  equations 
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(16) 


F"  +  -F'  +  £F  =  0 

r 


and 

G'  +  aX2G=  0,  (17) 

where  the  eigenvalue  k  remains  to  be  determined.  The  solutions  to  equations  16  and 
17  are 

F  =  AJ0(Ar)  +  BY0(Ar)  (18) 

and 

G  =  Ce-*at,  (19) 

where  Jo  is  the  Bessel  function  of  the  first  kind  of  order  zero,  Yo  is  the  Bessel  function 
of  the  second  kind  of  order  zero,  and  A,  B,  and  C  are  unknown  constants. 
Substituting  equations  18  and  19  into  equation  15  and  combining  constants  gives  the 
general  solution 

©  =  [EJ0{Pr)  +  DY0{Ar)]ellat ,  (20) 

where  E  =  AC  and  D  =  BC.  Applying  of  the  boundary  condition  at  r  =  0,  given  by 
equation  13,  yields 

D  =  0 .  (21) 

The  value  of  X  is  determined  by  substituting  equation  20  into  the  convective 
boundary  condition  (equation  14),  which  gives  the  following  expression 

=  BiJ0(k„R )  n  =  1,2,3,--,  (22) 

where  the  subscript  n  denotes  the  fact  that  an  infinite  number  of  eigenvalues  A.„ 
satisfy  equation  22,  J\  is  the  Bessel  function  of  the  first  kind  of  first  order,  and  the 
Biot  number  is  defined  as 


Bi  = 


hR 

k 


(23) 


Since  R  and  the  Biot  number  are  known  values,  the  eigenvalues  can  be  obtained  by 
solving  equation  22  numerically  or  graphically.  A  graphical  technique  was  used  and 
is  discussed  in  section  2.2.  With  the  eigenvalues  known,  the  general  solution 
becomes 


71=1 

An  expression  for  the  last  unknown  constant,  E„,  is  obtained  by  using  the  initial 
condition  (equation  12).  Substituting  equation  24  into  equation  12  gives 
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(25) 


©o  =2XA(V)- 

«  =  1 


An  explicit  expression  for  E„  can  be  obtained  by  multiplying  both  sides  of 
equation  25  by  rJo(X„r),  integrating  both  sides  from  r  =  0  to  r  =  R,  and  using  the  fact 
that  the  family  of  Bessel  functions  Jo(X„r)  are  orthogonal  with  respect  to  the 
weighting  factor  r  (see  Poulikakos  [1994]  for  details).  The  resulting  expression  is 


E.  =  - 


2©n 


D-  %R 

Bi  +  — — 


2\ 


Bi 


(26) 


Substituting  the  expression  for  E„  in  equation  26  into  the  general  solution  in 
equation  24  yields  the  final  equation  for  temperature: 


©  =  r-7;=2£z0o£ 

n=\ 


(27) 


2.2  Solution  Technique 

Solution  of  equation  27  for  penetrator  temperature  requires  a  defined 
one-dimensional  penetrator  geometry,  necessary  boundary  and  initial  conditions, 
appropriate  material  properties,  and  the  eigenvalues  that  satisfy  equation  22.  The 
geometry  considered  is  approximately  the  same  as  the  main  body  portion  of  the 
penetrator  used  in  Guidos  and  Weinacht  (1993),  depicted  in  Figure  2.  The  radius,  R, 
was  12.5  mm.  The  convective  heat  transfer  coefficient,  h,  adiabatic  wall  temperature 
temperature,  Tx,  and  initial  temperature,  T0,  were  1,700  W/m^K,  1285  K,  and 

295  K,  respectively  for  a  penetrator  velocity  of  approximately  1,500  m/s  (Guidos  and 
Weinacht  1993).  The  material  properties  used  for  each  material  are  listed  in  Table  1. 

Table  1.  Material  properties  for  WHA  and  DU  penetrators  (Weast  1975). 


k 

■ran 

c 

(W/m-K) 

mmmm 

(J/kg-K) 

WHA 

178 

17,650 

134 

DU 

25 

18,600 

117 

The  eigenvalues  needed  to  solve  the  infinite  series  (equation  27)  were  obtained  using 
a  graphical  technique.  This  process  involved  plotting  equation  22  with  X  as  the 
independent  variable  and  identifying  the  roots  of  the  equation  (i.e.,  recording  the 
values  of  X  where  equation  22  intersects  the  X  axis).  Figures  3  and  4  are  plots  of 
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Figure  3.  Graphical  depiction  of  characteristic  equation  for  extracting  eigenvalues  for  WHA. 


Figure  4.  Graphical  depiction  of  characteristic  equation  for  extracting  eigenvalues 
for  DU. 
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equation  22  for  WHA  and  DU,  respectively.  Since  the  solution  for  penetrator 
temperature  involves  an  infinite  series,  it  was  not  known  a  priori  how  many 
eigenvalues  would  be  needed  for  solution  convergence.  As  such,  the  first  20 
eigenvalues  for  each  material  type  were  identified  and  are  listed  in  Table  2.  As  will 
be  mentioned  later,  solution  convergence  occurred  rapidly  and  no  additional 
eigenvalues  were  needed. 


Table  2.  First  20  eigenvalues  satisfying  equation  22  for  WHA  and  DU  penetrators. 


WHA 

DU  I 

II 

^n 

^n  I 

1 

38.51 

94.21 

2 

309.02 

323.55 

3 

562.61 

570.81 

4 

814.82 

820.52 

5 

1066.61 

1070.98 

6 

1318.23 

1321.77 

7 

1569.76 

1572.73 

8 

1821.23 

1823.79 

9 

2072.66 

2074.92 

10 

2324.08 

2326.08 

11 

2575.47 

2577.28 

12 

2826.85 

2828.51 

13 

3078.23 

3079.75  j 

14 

3329.59 

3331.00  | 

15 

3580.96 

3582.26 

16 

3832.32 

3833.54 

17 

4083.67 

4084.82  | 

18 

4335.02 

4336.10 

19 

4586.37 

4587.39 

1  20 

4837.71 

4838.68 

With  the  eigenvalues  identified,  solving  equation  27  for  penetrator  temperature  is 
straightforward.  The  results  presented  in  the  next  section  were  obtained  using  the 
first  20  terms  of  the  infinite  series  expression.  Convergence  of  the  solution  was 
observed  to  occur  rapidly.  Calculated  temperature  was  found  to  change  less  than 
0.1  K  after  the  first  eight  terms  of  the  infinite  series  for  both  material  types. 


3.  Results  and  Discussion 


The  temperature  profile  across  the  penetrator  radius  was  calculated  with 
equation  27  for  each  material  at  several  flight  times.  The  profiles  are  plotted  in 
Figures  5  and  6  for  WHA  and  DU,  respectively.  The  temperature  increases  shown  in 


Temperature  Increase  (K)  3  Temperature  Increase  (K) 


lire  5.  Temperature  increase  inside  a  WHA  penetrator  for  various  flight  times. 


Figure  6.  Temperature  increase  inside  a  DU  penetrator  for  various  flight  times. 


both  plots  are  relative  to  ambient  (295  K).  The  temperature  increase  on  the  surface 
of  the  WHA  penetrator  reaches  a  value  of  300  K  after  3  s  of  flight.  With  a  flight 
speed  of  1,500  m/s,  this  corresponds  to  an  engagement  distance  of  4.5  km. 
Penetrator  flights  as  short  as  1  s  (1.5  km  distance)  result  in  surface  temperatures  of 
approximately  125  K  above  ambient.  An  even  shorter  flight  time  of  0.06  s,  which 
corresponds  to  a  typical  flight  distance  for  short-range  terminal  ballistic 
experiments,  results  in  a  temperature  increase  of  about  25  K.  The  relatively  high 
thermal  conductivity  of  WHA  results  in  very  little  thermal  variation  through  the 
radius,  thus  the  entire  body  of  the  penetrator  is  heated  to  elevated  temperatures. 

The  surface  temperature  increase  for  penetrators  made  of  DU  were  higher  than 
those  calculated  for  WHA.  The  temperature  increase  in  DU  after  3  s  was  400  K, 
compared  to  300  K  for  WHA.  The  interior  region,  however,  was  cooler  than  the 
corresponding  WHA  case.  This  is  due  to  the  low  thermal  conductivity  of  DU,  which 
is  lower  than  WHA  by  a  factor  of  approximately  7.  As  a  result,  there  exists  a 
significant  thermal  gradient  in  the  DU  penetrator.  This  was  evident  even  for  the 
shortest  flight  time  considered,  0.06  s,  which  resulted  in  a  temperature  increase  of 
55  K  on  the  surface  and  no  temperature  rise  2.5  mm  below  the  surface. 

The  ductile-brittle  transition  temperature  of  WHA  is  approximately  room 
temperature  under  quasi-static  loading  conditions  (Larson  and  DeLai  1980;  Gurwell 
1986).  The  transition  temperature  will  increase  for  high-rate  loadings  (Freund  1990), 
but  it  is  currently  unknown  how  much  the  increase  will  be.  As  such,  the  calculated 
temperatures  for  the  WHA  penetrator  shown  in  Figure  5  suggest  the  mode  of 
fracture  during  a  terminal  ballistic  event  may  vary,  depending  on  the  flight  time  or 
distance.  Short-range  engagements  may  result  in  the  penetrator  undergoing  brittle 
failure,  whereas  longer  flight  times  may  produce  ductile  failure.  Attempts  to 
characterize  and  model  the  interaction  of  the  penetrator  with  a  complex  armor  that 
applies  lateral  loading  will  likely  need  to  account  for  the  difference  in  fracture 
behavior  of  the  penetrator  material. 

The  temperature  distribution  in  the  DU  penetrator  offers  a  different  challenge  to 
modeling  its  terminal  ballistic  performance.  For  this  material,  a  significant  thermal 
gradient  develops  through  the  radius  of  the  penetrator.  Although  the  mode  of 
fracture  will  likely  remain  ductile  because  of  the  relatively  low  transition 
temperature  (Eckelmeyer  1991),  a  crack  originating  at  the  surface  of  the  penetrator 
will  propagate  into  material  that  has  a  decreasingly  lower  yield  stress.  Berger  et  al. 
(1993)  performed  crack  propagation  experiments  with  steel  that  had  thermal 
gradients  through  the  thickness,  where  the  crack  propagated  into  higher 
temperature  material.  The  experiments  showed  fracture  toughness  to  increase  with 
material  temperature,  thereby  leading  to  crack  arrest  when  it  propagated  to  material 
with  a  high  enough  temperature.  Similar  effects  will  likely  be  present  for  DU 
penetrators  at  intermediate  and  long  flight  times,  except  in  this  case  the  gradient  will 
promote  unstable  crack  growth. 
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4.  Conclusions 


A  one-dimensional  transient  heat  conduction  model  was  derived  to  provide  the 
temperature  distribution  inside  a  KE  penetrator  made  from  either  WHA  or  DU 
traveling  through  air  at  1,500  m/s.  The  governing  equation  was  solved  using  the 
separation  of  variables  technique,  resulting  in  an  infinite  series  expression  for  the 
penetrator  radial  temperature  distribution.  A  flight  time  of  3  s  was  sufficient  to  raise 
the  temperature  of  the  penetrator  surface  by  300  and  400  K  for  WHA  and  DU, 
respectively.  The  surface  temperature  increase  calculated  for  a  much  shorter  flight 
time  of  0.06  s  was  25  K  for  WHA  and  60  K  for  DU.  The  high  thermal  conductivity 
for  WHA  resulted  in  a  somewhat  uniform  radial  temperature  distribution.  The  heat 
absorbed  at  the  penetrator  surface  was  quickly  conducted  to  the  interior  of  the 
penetrator  body.  In  contrast,  the  low  thermal  conductivity  of  DU  suppressed 
thermal  conduction,  resulting  in  a  high  surface  temperature  with  a  significant 
thermal  gradient  across  the  penetrator  radius. 

The  temperatures  calculated  in  this  study  are  significant  from  the  perspective  of 
modeling  the  terminal  ballistic  interaction  of  WHA  and  DU  penetrators  with 
complex  armors,  particularly  when  fracture  is  considered.  The  proximity  to  the 
ductile-brittle  transition  temperature  to  the  temperatures  calculated  for  the  WHA 
case  indicate  that  fracture  mode  transition  must  be  considered  when  developing 
failure  models.  Additionally,  die  thermal  gradients  calculated  for  the  DU  penetrator 
may  have  significant  implications  regarding  dynamic  crack  propagation,  which  may 
also  need  to  be  accounted  for  in  detailed  failure  models. 
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